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We investigate the orbital diamagnetism of a weak-doped bilayer graphene (BLG) in spatially 
smoothly varying magnetic field and obtain the general analytic expression of the orbital suscepti- 
bility of BLG, with finite wave number and Fermi energy, at zero temperature. We find that the 
magnetic field screening factor of BLG is dependent with the wave number, which results in a more 
complicated screening behavior compared with that of monolayer graphene (MLG). We also study 
the induced magnetization, electric current in BLG, under nonuniform magnetic field, and find 
that they are qualitatively different from that in MLG and two-dimensional electron gas (2DEG). 
However, similar to the MLG, the magnetic object placed above BLG is repelled by a diamagnetic 
force from BLG, approximately equivalent to a force produced by its mirror image on the other side 
of BLG with a reduced amplitude dependent with the typical length of the systems. BLG shows 
crossover behaviors in the responses to the external magnetic field as the intermediate between MLG 
and 2DEG. 

PACS numbers: 75.20.-g, 73.63.-b, 75.70.Cn, 81.05.Uw 



1. INTRODUCTION 

Bilayer graphene (BLG), as a significant graphene-related material, has attracted much attentions [l--)] due to 
its unusual electronic structure. Formed by stacking two monolayer graphene (MLG) in Bernal stacking, Bilayer 
graphene has four inequivalent sites in each unit cell, including Ai and Bi atoms on the top layer and A2 and B2 
atoms on the bottom layer. The distance a between Ai and Bi is about 0.142 nm and the vertical separation of the 
two layers d is about 0.334 nm. On its band structure side, around the point where the conduction and valence band 
touch, BLG has a quadratic energy dispersion [1] similar to the regular two-dimensional electron gas (2DEG) but its 
low-energy effective Hamiltonian is chiral without bandgap similar to the MLG. Another unique feature of BLG is 
that a widely tunable bandgap can realize conveniently by introducing an electrostatic potential bias between the up 
and bottom layer [1, 4-8]. 

The magnetic susceptibility of electronic systems comes from two contribution: One is Pauli paramagnetism which 
stems from spin polarization; The other is Landau diamagnetism which stems from the circulation of orbital currents. 
The orbital diamagnetism of carbon systems has attracted the interest of both experimental and theoretical physicists 
for a long time. It was firstly found by Krishnan [!)] that the diamagnetic susceptibility of the bulk graphite is large and 
anisotropic. McClure [10] showed it arises from the Landau quantization of 2D massless Dirac fermions, which results 
in a delta function peak at zero energy in the orbital diamagnetism of graphite. After the experimentally fabrication 
of graphene, a great deal of works have concerned about the orbital magnetism of graphene-related systems, such 
as nodal fermions [11], disordered graphene [12-14], few-layered graphene [15-17], graphene in nonuniform magnetic 
field [18, 19]. Quite recently, novel paramagnetic susceptibility has been found in MLG and BLG with bandgap [20] 
and doped graphene to the first order in the Coulomb interaction [21]. 

Safran [22], Koshino and Ando [15] have derived the analytical expression of orbital susceptibility x(0) for bilayer 
graphene, but their results are just limited to the case of zero wave number. In this article we generalize them to 
general cases and analytically study the orbital diamagnetism of a weak-doped BLG (i.e., the Fermi energy e^? is much 
smaller than the interlayer hopping energy.) in nonuniform magnetic fields. 

This article is organized as follows. In Sec. II, the effective Hamiltonian of bilayer graphene and its corresponding 
eigenstates and eigenenergies are introduced. In Sec. Ill, the orbital susceptibility of bilayer graphene is studied. In 
Sec. IV, we investigate and discuss the responses of bilayer graphene to several specific external magnetic field. The 
conclusion is given in Sec. V. 
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2. BILAYER GRAPHENE EFFECTIVE MODEL 



111 the low energy and long wave regime, the BLG Hamiltoiiian near a K point, in the absence of a magnetic field, 
can written as an excellent approximate form [1] (we set c = 1 = in this paper): 
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where 7 = 3ta/2 « IQ^m/s is the monolayer graphcnc Fermi velocity, t k. 2ieV is the in-plane hopping energy. 



a « 0.142 nm is the in-plane interatomic distance, A « Q.ibeV is the interlayer hopping energy, k = {k^, ky) 



-iV 



is a 2D wave-vector operator, k± = k^ iky. In order to solve the eigen equation of the Hamiltonian (1), the wave 
function can be expressed as (V'An V'Sn V'A2j "032)1 where the four components represent the Bloch functions at Ai, 
Bi, A2 and B2 sites, respectively. Follow the stipulation of Ando ["J-!], and define 



6(fc) 



jk ~ e{k)smip, 
~ e{k)cos^l>. 



Then the corresponding cigenstates of Eq.(l) arc given as 
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where is the area of the system, s = -1-1 and —1 denote the conduction and valence bands, respectively, j = 1 
and 2 specifies two subbands within the conduction or valence bands, 6k — arctan(fcj,/fc2;) is the polar angle of the 
momentum fc, 
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and 



/ scos(V'/2) \ 



sin(VV2) 
-s sin(V'/2) 

V - cos(v/2) ; 

The corresponding eigenenergies of Eq.(l) are 
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V sin(^/2) ; 



(7) 



2se(fc)sin2(?A/2), 
2se(fc)cos2(V'/2). 
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Considering a magnetic field B{r) = [V x A{r)]z, the Hamiltonian for tlie system is: H = Hq + Hi, with 
Hi ~ ~ J cPrja{r)Aa{r,t). The current operator at Tq is given 

6 

jairo) = ^[vaSir - ra) + 5{r - ro)va], {a = x,y) (9) 

where Va is velocity operator 

(a = .T,y) (10) 
ax,y are the Pauli matrices which act on the sublattice space within a layer. 
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3. ORBITAL SUSCEPTIBILITY 



The finite wave number susceptibility x{q) can be obtained through the Kubo formula [24]. Within the linear 
response theory, the external vector potential A and the its induced 2D electric current density j have a relation 



j,.{q)=^K,,,{q)A,{q), 



(11) 



and the orbital susceptibility x{q) and the response tensor K^^{q) are related by 



-)• 



In the first order perturbation, we have 



ss'jj'k J ^ ' J \ / 



(12) 



(13) 



where g = g^gs = 4 is the total degeneracy, k' = k + q, esj{k) is the eigenenergy given by Eq.(8), /(e) is the 
Fermi distribution function /(e) = [1 + cxp/3(e — e_F)]~-^ where ep is the Fermi energy, j3 = l/{kBT). Iss'jj' is the 
current-current response matrix element expressed by 



(14) 



which determines the weight of contribution of the transition from subband j to j' , with ss' = +1 and —1 denote the 
intraband and interband transition, respectively. 

Define the effective mass m = A/(27^) sa 0.033me and the Fermi wave number kp = ^2mtF- For a low-energy 
theory of bilayer graphene, there is a natural high-energy cutoff wave number A = V2mA = A/7. When /cp, g <C A, 
i.e., when bilayer graphene is weak-doped and the external field is smooth enough compared to the cutoff wavelength, 
the transitions intra the same j subband dominate the contribution to the response function. One of these transitions 
is the transition intra the j — \ subband: 



9 ^ / [e.i(fc)]-/[e,a(fcO] 
e,i(fc)-e,.i(fc') 
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(15) 



In the low energy limit, it can be approximately given as 



with ek ~ k'^ / {2m) and 



ss'h 



(16) 



Ff,y = {Sf^iS^i - Sf,2Su2) [fc^cos26ifc/ + k'^cos2ek + 2fcfc'cos(6ifc -t^ Ok')] 
+ {S^,i6^2 + S^,2S^i) [k^ sin2ek' + k'^ sm2ek + 2kk' smiOk + ek')] . 



(17) 



4 



The other is the transition intra the subband j = 2: 



if..,22(g) --Jll^_ ..2(fc)-6.s'2(fc') 



ss'k. 



(18) 



By using the fact that the subband j = 2 in the conduetion band (s = +1) is empty in the weak-doped limit, it ean 
be approximately given as 



a: 



2ge^ 
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(19) 



with 



G^,^ = {5^,i5^i - 5^25u2) [fc^ cos26'fc + fc'^ cos26'fc/ + 2kk' cos{0k + Ok')] 

+ {S,,iS,2 + 5^,25ui) [fc' sin20fc + k'^ sin20fc, + 2kk' sm{0k + Ok')] . (20) 
From the above, we firstly obtain the analytic expression of the susceptibility of bilayer graphene at zero temperature 



as: 
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(21) 



where 9{x) is the step function defined by 9{x) = l{x > 0) and 0(a; < 0), which is a centra result of our work. From 
Eq.(21), on the one hand, we can give the orbital susceptibility of BLG at zero wave number as 



Xiq = 0;ei?) 



Bttto 



A + 3 



(22) 



This result is same as the result given in [15, 22], which shows a logarithmically diverging behavior at ep = 0. On 
the other hand, we can give the orbital susceptibility of BLG at zero Fermi energy as 



x(q;eF = 0) 



loa 



2A 



(23) 



This also shows a logarithmically diverging behavior at g = 0. It is easy to see that just by replacing ep with 7(7/2 
and increasing the coefficient to its double times, we can transform the susceptibility x(g = 0; e^?) into xiQi = 0). 

The orbital magnetic susceptibility x{q) of bilayer graphene as a function of wave number is shown in Fig.[l]. For 
comparing the wave number-dependent behaviors of susceptibility of the MLG, BLG and 2DEG systems, we provide 
below the finite wave number susceptibility of MLG and 2DEG, which has been given in Ref.[lS], 
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x(q;eF) = 
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(25) 



At g = 0, the susceptibility of BLG — x(0, ei?) oc log(A/ei?) is rather different from that of MLG, where ~x{^i ^f) oc 
6{€f), and 2DEG, where — x(0,ei?) oc 1. By comparing their diverging behaviors, it can be found that the BLG in 
some sense shows an intermediate behavior between the MLG and 2DEG. For small q, the — x(g; <^f) of BLG deviates 
from the — x(0;ei?) as {q/2kF)^, and falls more rapidly as q increase. On the other hand, the susceptibility of MLG 
vanishes while 2DEG maintains as a constant for the whole regime q < 2kF (see Fig.[l] of Ref. [18]). At g = 2kF, the 
susceptibility x(2fci?; ep) of MLG and 2DEG are both constants (zero for MLG) which are independent of the Fermi 
wave number kp, but for BLG, we have 



X{2kp;tp) 



ge" 
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which is dependent of the kp. In contrast to the MLG, the susceptibility of BLG has no singular behavior at q ~ 2kp, 
and it is continuous as well as its first derivative. For large g, especially for q ^ 2fci?, —xi.l) of BLG rapidly approaches 
the curve of Eq. (23) and falls as \og{l/q), very different from that of MLG where the susceptibility falls off more 
rapidly l/g) and 2DEG where it falls as l/g^. Due to having the same parabolic energy dispersion in the low 
energy limit, susceptibility of Bilayer graphene and 2DEG share the same term ge^ {1 — Akj;' / q^)'^^^6{q — 2kF) / {^ATmi) . 



4. RESPONSES TO SPECIFIC EXTERNAL MAGNETIC FIELD 



Now we study the responses of bilayer graphene to different types of magnetic field. First let us consider the case of 
a neutral BLG (i.e., ep = 0) under a sinusoidal magnetic field B{r) = Bq cosgxe^. Defining x{<l) = xiQj = 0), we 
have the induced magnetization m(r) = x(g)B(r), the induced current j{r) = qx{q)BQ s'mqxey, and the z component 
of induced counter magnetic field on BLG 



Bindir) = -ag{q)B{r), 



with the magnetic field screening factor 



ag{q) 



2m 



log 



2A 



(27) 



(28) 



When q — > 0, we have ag{q) — > 0; i.e., under a constant magnetic field there is no counter magnetic field on BLG, 
which is different from that of MLG. In MLG, the magnetic field screening factor is fixed and independent of q and 
the specific form of external field, the induced magnetic field above the graphene layer is simply equivalent to the field 
of a mirror image of the original object reflected with respect to the graphene layer but reduced by Ug [18]. However, 
this argument does not fit for the case of BLG, since the magnetic field screening factor of BLG is dependent of q and 
complicated. 

Next we consider the case of a line current / fiowing along the +y direction above the BLG, and passing through 
the point (0,0,6?) {d > 0). The z component of the magnetic field on BLG is B{r) = —2Ix/{x^ + dP). For ep = 0, 
the induced magnetization can be given as 



m(r) = — 



Ige^ 



27rm + dp 



log(2A-\/ x'^ + d"^) + — arctan ^ — 7b + ^ 
X d 3 



(29) 



here7£; « 0.577 is the Euler constant. Byusing ji„d = Vxm(r), we obtain the induced electric current ji„d('") = jy^y, 
where 



3v 



1 



27rm (x2 + rf2)2 



{d^ - x^) 



log(2Av/a;2 +d2) - + 1 



- 2dx arctan - + {x^ + d^) 
d 



(30) 
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The integral of jy in x exactly equals to 0, which means the external electric current I can not induce an effective 
transport electric current on the BLG, in contrast to that of MLG, where / induces an effective electric current —ctgl- 
The induced magnetic field on BLG can be given as B^ndif) ~ B^e^ with 



s jx'^ + (P- _ 1 



2dx I log ' " + 7_E - ^ 1 + (d^ - x^) arctan ^ 



(31) 



In MLG, the induced magnetic field —aglx/{x'^ + (P) is equivalent to the field created by a current —agl flowing at 
z = —d. This argument does not fit for the case of BLG as shown in Eq.(31). At large distance x ^ d, the induced 
magnetic field is proportion to ~ 1/x'^ comparing with ^ 1/a; in MLG. However, the original current is repelled by a 
force « ag{l/2d)I^ /d per unit length, which can be approximately but not exactly considered as a force created by 
a current ag{q)I at z = —d with q = l/2d. 

As another typical example, we study the the induced magnetization, electric current and magnetic field by a 
magnetic monopole Qm laying above the BLG. Suppose Qm is located at the point (0, 0,(i), {d > 0), and the BLG 
plane is z = 0. The magnetic field perpendicular to the BLG is given as B{r) = q^d/ir"^ with r — \J x"^ -\- y"^. 

For neutral BLG, the induced magnetization is given by 



-3/2 



where the function 

= — / zJo(z)logze-^/^dz. (33) 
Jo 

At small distance (r <C d) and large distance (r ^ d), the induced magnetization can be written as 



{log2Ad + 7B-|-f(r/rf)2 [log2Ar-i]}/d2 (for r < d) 
l/r^ (for r > d) 



m{r) = X { (34) 



(35) 



It is interesting to see that at large distance the induced magnetization of BLG m(r) cx l/r^, while that of MLG 
cx l/r and 2DEG cx l/r^; i.e., the BLG shows a behavior as the crossover from MLG to 2DEG. The integral of m(r) 
over the plane has a logarithmically diverging behavior cx logi? when the distance i? — > cx3. According to Eq.(32), 
the corresponding electric current can be given by j{r) ~ -~{drn/dr)eg = jgeg. When r is small or large, we obtain 

2 r 3rlog2Ar/d4 {for r « d) 
3e = -; X <^ (36) 



ATTm 



(for r > d) 

(37) 



Remind that our results are confined to the limit Ar 3> 1, and therefore jr is positive through out the realistic 
distance. The current jg in MLG oc r/(r^ + d'^)^^, whose asymptotic form is cx r at small distance and cx l/r^ at 
large distance. We can find that the induced currents jg of MLG and BLG are qualitatively different in all distance. 
However, similar to the case of line current, the force between the monopole and the BLG can be approximately 
written as ag{q)q^/{2dY with q — 1/d, which has the same form as that of MLG agq^/{2d)'^ . 

For doped bilayer graphenc [ep ^ 0), the corresponding response current jg{r) for different values of fc^^ are shown 
in Fig. [2]. We can find that the current change slightly as the Fermi wave number increase from kpd = Q to kpd = 1] 
i.e., the induced current has a weak Fermi wave number dependence, which is rather different from that of MLG. This 
arises from the fact that the dominant term involving kp in the susceptibility is logarithmically. 

Different from the traditional 2DEG, MLG has a peculiar property [IS]: (1) The counter field induced by the 
response current mimics a mirror image of the original object. (2) The object is repelled by a diamagnetic force from 
the MLG, as if there exists its mirror image with a reduced amplitude on the other side of MLG. With the investigation 
above, we find that argument (1) can be not extended to BLG. However, the argument (2) still be approximately 
correct for weak-doped BLG, and it only needs to replace the constant reduced amplitude with a reduced amplitude 
dependent with the typical length of the systems. The BLG, in some sense, still show an intermediate behavior 
between the MLG and 2DEG. 



7 



4 




-1' ■ ■ ■ ■ ' 

1 2 3 4 5 

r/d 

FIG. 2: Electric current je(r) on bilayer graphene induced by a magnetic charge qm a.t z = d. Here we use A/kp = 10. 

It is significant to compare the contribution of Landau diamagnetism to the whole magnetism with that of Pauli 
paramagnetism in BLG. At = 0, the Pauh spin susceptibility which is given by the density-density 

response function [25], is equivalent to g„m^^ log4/27r, where is the Bohr magneton. With Eq.(23), we obtain 
the ratio ^^pm/^orb ^ (0.03)V(log 2A/g), which is rather small in our theory {for q ^ A). 

The temperature also has an influence on the BLG diamagnetism. For q ~ 0, this has been discussed by Safran 
[22], who shows the susceptibility with finite temperature and Fermi energy can approximately take the form 

x(9 = 0;ej.;r)cxlog|Ai-/A|-l + (/i+/fcBr)log|AH/M-l (for « A). (38) 

here /i± = e^? ± kBT/2. For finite q, we expect that xil) deviates from logq in regime q < ^1m\tF — fc_Br/2|, which 
means the temperature has notable affection on the diamagnetism of a neutral BLG when the typical length scale 
exceeds 2tt / y/2mkBT , about TOfim at T = IK. 

All of the above, we assume an ideal 2D BLG electron gas and ignored the distance d between up layer and bottom 
layer of bilayer graphene. In order to obtain an analytic expression of the susceptibility of bilayer graphene, we assume 
the wave number kp and q is much smaller than the cutoff wave number A. An alternative farther investigation can 
be focused on a more general case that fc^^ is commensurate to or even larger than A. 

5. CONCLUSIONS 

In this article, we study analytically the orbital magnetic susceptibility of a weak-doped bilayer graphene (BLG) in 
spatially smoothly varying magnetic fields by the low-energy Hamiltonian and obtain the general analytic expression 
of the orbital susceptibility of BLG, with finite wave number and Fermi energy, at zero temperature. The induced 
magnetization, electric current by the nonuniform magnetic fields in BLG are studied which are different from that 
of MLG and 2DEG, but the argument, that the magnetic object placed above the BLG is repelled by a diamagnetic 
force which is equivalent to a force produced by mirror image on the other side of BLG, still be approximately hold, 
only by replacing the constant reduced amplitude with a reduced amplitude dependent with the typical length of 
the system. Logarithmically-dependent behaviors are found extensively exist in both the orbital magnetism and the 
induced physical quantities by specific external field. BLG shows crossover behaviors in the responses to the external 
magnetic field as the intermediate between MLG and 2DEG. The weak Fermi wave number dependent behaviors, as 
a distinctive electric property of BLG, arc found in induced magnetization and electric current. 
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